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Abstract 

The kinetic energy of barotropic flow coupled to an infnitely massive ro- 
tating sphere by an unresolved complex torque mechanism is approximated 
by a discrete spin-lattice model of fluid vorticity on a rotating sphere, analo- 
gous to a one-step renormalized Ising model on a sphere with global interac- 
tions. The constrained energy functional is a function of spin-spin coupling 
and spin coupling with the rotation of the sphere. A mean field approxi- 
mation similar to the Curie- Weiss theory, modeled after that used by Bragg 
and Williams to treat a two dimensional Ising model of ferromagnetism, is 
used to find the barotropic vorticity states at thermal equilibrium for given 
temperature and rotational frequency of the sphere. A fixed point equation 
for the most probable barotropic flow state is one of the main results. 

This provides a crude model of super and sub-rotating planetary atmo- 
spheres in which the barotropic flow can be considered to be the vertically 
averaged rotating stratified atmosphere and where a key order parameter is 
the changeable amount of angular momentum in the barotropic fluid. Using 
the crudest two domains partition of the resulting fixed point equation, we 
find that in positive temperatures associated with low energy flows, for fixed 
planetary spin larger than Qc > there is a continuous transition from a 
disordered state in higher temperatures to a counter-rotating solid-body flow 



state in lower positive temperatures. The most probable state is a weakly 
counter-rotating mixed state for all positive temperatures when planetary 
spin is smaller than Qc- 

For sufficiently large spins fl > 2flc, there is a single smooth change from 
slightly pro-rotating mixed states to a strongly pro-rotating ordered state as 
the negative value of T increases (or decreases in absolute value). But for 
smaller spins Q < 2Qc there is a transition from a predominantly mixed state 
(for T ^ — 1) to a pro-rotating state at l3^{Q) < plus a second /?^^\ for 
which the fixed point equation has three fixed points when {3 < {3^ instead 
of just the pro-rotating one when < [3. An argument based on comparing 
free energy show that the pro-rotating state is preferred when there are three 
fixed points because it has the highest free energy - at negative temperatures 
the thermodynamically stable state is the one with the maximum free energy. 

In the non-rotating case {VL = 0) the most probable state changes from a 
mixed state for all positive and large absolute-valued negative temperatures 
to an ordered state of sohd-body flow at small absolute-valued negative tem- 
peratures through a standard symmetry-breaking second order phase tran- 
sition. The predictions of this model for the non-rotating problem and the 
rotating problem agrees with the predictions of the simple mean field model 
and the spherical model. 

This model differs from previous mean field theories for quasi-2D turbu- 
lence in not fixing angular momentum nor relative enstrophy - a property 
which increases its applicability to coupled fluid-sphere systems and by ex- 
tension to 2d turbulent flows in complex domains such as no-slip square 
boundaries where only the total circulation are fixed - as opposed to classical 
statistical equilibrium models such as the vortex gas model and Miller- Robert 
theories that fix all the vorticity moments. Furthermore, this Bragg mean 
field theory is well-defined for all positive and negative temperatures unlike 
the classical energy-enstrophy theories. 
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1 Introduction 



Previous statistical equilibrium studies of inviscid flows on a rotating sphere 
with trivial topography focussed on the classical energy-enstrophy formula- 
tion [17J of PDE models such as in [16] did not find any phase transitions [H] 
to super-rotating solid-body flows. The main reason is that these Hamilto- 
nian PDE models which conserve the angular momentum of the fluid, are not 
suited to the study of a phenomena like the super-rotation of the Venusian 
middle atmosphere that depends on transfer of angular momentum between 
the atmosphere and the solid planet. 

We therefore consider instead a single layer of barotropic fluid of fixed 
height coupled by a complex torque mechanism to a massive rotating sphere 
as a simplified model of planetary atmosphere. The approach in this paper is 
to consider different possible thermodynamic properties of the steady states 
of coupled planetary atmospheres as mean field statistical equilibria - through 
a natural relation of fluid vorticity systems to electromagnetic models - with 
the aim of showing explicitly that sub-rotating and super-rotating flow states 
can arise spontaneously from random initial states. 

A simple mean field method on a related model [19], (where spins are al- 
lowed to take on a continuous range of values), has successfully found phase 
transitions for barotropic flows coupled to a rotating sphere. There as here, 
angular momentum of the fluid is allowed to change; and the expected val- 
ues of total circulation and relative enstrophy are fixed in this simple mean 
field theory. This leads to the exciting result that sub-rotational and super- 
rotational macrostates are preferred in certain thermodynamic regimes. The 
Bragg mean field theory presented in this paper corroborates some of these 
previous results but has the important new property that relative enstrophy 
is constrained only by an inequality while total circulation is held fixed at 
zero. Further discussion of these important points will follow in the next 
section. 

Our results are also in agreement with Monte Carlo simulations on the 
logarithmic spherical model in the non-rotating case [T^] and the rotating 
sphere [T], The phrase spherical model refers to the fact that under a 
discrete approximation, the microcanonical constraint on the relative enstro- 
phy becomes a spherical constraint first introduced by Kac [18]. We shift 
gears here from the spherical model used by Lim - for which closed form so- 
lutions of the partition function required considerable analytical effort [12], 
[S] - to a simplified discrete spin model which enables us to use the Bragg 
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method to approximate the free energy, which in turn allows us to find an 
interesting analytical solution to the coarse grained stream function for the 
resulting mean field theory. While the overall qualitative agreement between 
this Bragg model, Lim's exact solution of the spherical model and the simple 
mean field model mentioned above is good, there are subtle differences in the 
detailed predictions of these models. They are however all within the scope 
of applications to super and sub-rotations in planetary atmospheres such 
as on Venus, chiefiy because current observations are not refined enough to 
distinguish between them. It will be necessary to perform detailed direct 
numerical simulations of coupled geophysical fiows in order to gather enough 
data for a future comparative study which can choose between these three 
models. Nonetheless, we will discuss in the next section how purely theo- 
retical considerations of the complex boundary conditions that exist in the 
coupled barotropic fiuid -sphere system and in a class of quasi-2d turbulent 
fiows in complex domains allow us to choose the current set of spin-lattice 
/ Eulerian vorticity models with minimal constraints over the vortex gas / 
Lagrangian models [13], [S] and Miller- Robert theories [S], [TU]. 

Our discrete vorticity model is a set X = {afj} of 2N random sites on 
the sphere with uniform distribution, each with a spin Si G {+1,-1}, and 
interaction energy with every other site as a function of distance. Each 
spin also interacts with the planetary spin which is analogous to an external 
magnetic field in the Ising model - the sum of this interaction is proportional 
to the varying net angular momentum of the fiuid relative to a frame that 
is rotating at the angular velocity of the solid sphere. Contribution to the 
kinetic energy from planetary spin varies zonally, however, and this makes the 
external field inhomogeneous and difficult to deal with analytically. Bragg 
and Williams [2] used a one step renormalization to investigate properties 
of order-disorder in the Ising Model of a ferromagnet. As the discretized 
model for barotropic fiows coupled to a rotating sphere is similar to the Ising 
Model of a ferromagnet in an inhomogeneous field, we use the Bragg- Williams 
renormalization technique to infer the order-disorder properties of the fiuid. 

It is a well established fact [TB] that systems of constrained vortices have 
both positive and negative temperatures, and the present model is no excep- 
tion. In this setup, using the simplest two domains partition of the surface of 
the sphere, we find a positive temperature continuous phase transition to the 
sub-rotating ordered state for decreasing temperatures if the planetary spin 
is large enough. There is no evidence of a phase transition for a non-rotating 
sphere in positive temperature. We also find transitions to a super-rotating 
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ordered state when the negative temperature has small absolute values - that 
is for very high energy. 

2 Statistical mechanics of macroscopic flows 

Here we summarize some of the pertinent points concerning the application 
of a statistical equilibrium approach to macroscopic flows. The largely 2-d 
flows we are concerned with in this paper are in actuality non-equilibrium 
phenomena even in the case of nearly inviscid flows - an assumption that is 
mainly correct for the interior of geophysical flows. The main reason that 
equilibrium statistical mechanics is applicable here - the existence of two 
separate time scales in the microdynamics of 2D vorticity - is best formulated 
as the well-known physical principle of selective decay which states that the 
slow time scale given by the overall decay rates of enstrophy and kinetic 
energy in damped unforced flows is sufficiently distinct from the fast time 
scale corresponding to the inverse cascade relaxation of kinetic energy from 
small to large spatial scales so that several relaxation periods fits in a unit 
of slow time. That means that the total kinetic energy (and enstrophy) may 
be considered fixed in the time it takes for the eddies to reach statistical 
equilibrium. Furthermore, the principle of selective decay states that the 
asymptotic properties of the damped 2d fiow is characterized by a minimal 
enstrophy to energy ratio which is dependent on the geometry of the fiow 
domain. One of the key properties of enstrophy known as the square-norm 
of the vorticity field then implies that these minimum enstrophy states are 
associated with large scale ordered structures such as domain scale vortices. 
In the case of the specific problem in this paper, these coherent structures 
are super- and sub-rotating solid-body fiow states. 

Aside from the above general points on the applicability of equilibrium 
statistical mechanics to 2D macroscopic fiows, the specific properties of a 
given fiow problem [6] will decide which one of the many statistical equi- 
librium models [13], [8], [17], [9], [lO] is suitable. The coupled nature of 
the fiow problem in this paper where a complex torque mechanism transfers 
angular momentum and energy between the fiuid and rotating solid sphere, 
clearly eliminates all Lagrangian or vortex gas models [13], [8] on the surface 
S"^ because they conserve angular momentum - a result that follows directly 
from Noether's theorem. Moreover, since none of the vorticity moments are 
conserved in the coupled fiows except for total circulation which is fixed at 



5 



zero by Stokes theorem, independently (or in spite) of the couphng between 
fluid and solid sphere, the Miller- Robert class of models P, [ID] cannot be 
used here because they are designed to fix all the vorticity moments. It is 
therefore not surprising that a significant discrepancy was found between 
statistical equilibrium predictions based on a vortex gas theory and exper- 
imental (and direct numerical simulations) of decaying quasi-2D flows in a 
rectangular box with no-slip boundary conditions [6j. 

What remains are the classical Kraichnan models - also known as absolute 
equilibrium models - and variants such as those used in ^16j. The problem 
with these are two-fold - (i) they model a 2d fluid flow that is uncoupled to 
any boundary or has periodic boundaries and hence angular momentum is 
held fixed - clearly unsuitable for our aims here, and (ii) they are equiva- 
lent to the ubiquitous Gaussian model and can be easily shown to not have 
well-defined partition functions at low temperatures. One solution taken up 
in [12], [5] is to impose instead of a canonical enstrophy constraint, a micro- 
canonical one which leads to a version of Kac's spherical model. Combined 
with zero total circulation and allowing angular momentum to fluctuate - 
that is coupling the barotropic fluid to an infinitely massive rotating solid 
sphere- this approach yielded exact partition functions and closed form ex- 
pressions for phase transitions to self-organized (or condensed) super- and 
sub-rotating flows. As expected, these analytic results agree perfectly with 
Monte-Carlo simulations of the spherical model in [TB], PQ, An argument 
similar in style to the above discussion of the Principle of Selective Decay 
was needed to justify the choice of a microcanonical (and hence spherical) 
constraint on relative enstrophy in this class of coupled flows while the zero 
total circulation constraint needs none since it is implied by topological ar- 
guments - Stokes theorem on the sphere. The spherical model unlike the 
Bragg model in this paper is not based on a mean field assumption. This can 
be a significant advantage only in modelling more complex geophysical flows 
where large deviations techniques cannot be used to establish the asymptotic 
exactness of the mean field. 

One of the aims of this paper is therefore to relax the enstrophy constraint 
- that is, we impose neither a canonical enstrophy constraint like in [T7] nor 
a microcanonical enstrophy constraint like in [12j, [15j, [Ij, [?j - and 
still derive a physically sound mean field theory for a comparatively simple 
geophysical problem. We will show that the relative enstrophy in the Bragg 
mean field theory is constrained instead by an upper bound. 

Another aim in this paper is to show that the Bragg method can be ex- 
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tended to the situation as in this coupled flow where the kinetic energy is 
just a functional and not a Hamiltonian of the barotropic flow. This is a 
subtle but relevant point because the original Bragg method formulated for 
the Ising model of ferromagnetism worked with a Hamiltonian. The first 
author reported in |5] and [1], [1] that the kinetic energy of the fluid compo- 
nent of the coupled fluid - rotating sphere system cannot be a Hamiltonian 
for the evolution of the vorticity field because otherwise its 5*0(2) symme- 
try - a property that is easily shown - implies the conservation of angular 
momentum of the fluid component. 

3 Coupled Barotropic Fluid - Rotating Sphere 
model 

We are interested in the physical possibility of planetary atmosphere at ther- 
mal equilibrium with the rotating solid planet - the massive sphere is treated 
as an infinite reservoir for both kinetic energy and angular momentum of the 
barotropic (vertically averaged) component of the fluid. Thus we consider 
the atmosphere as a single layer of incompressible inviscid fluid interacting 
through an unresolved complex torque mechanism with the solid sphere. This 
is equivalent to formulating a Gibbs canonical ensemble in the energy of the 
fluid with an implicit canonical constraint on the angular momentum of the 
fluid. More realistic models include effects of forcing and damping, from solar 
and 'ground' interaction as well as internal friction. 

As our interest is in global weather patterns, we will consider the problem 
on the spherical geometry, with non-negative planetary angular velocity Q > 
0, and relative stream function ip - a characteristic of 2-d flows; q describes 
the total vorticity on the sphere, 

q = uj + 2ncos (6) 

where the first term u is the vorticity relative to the rotating frame in which 
the sphere is at rest, and the second gives planetary vorticity due to the 
angular velocity of the sphere in terms of the rotational frequency Q and the 
co-latitude 9. We note the fact that a distinguished rotating frame with fixed 
rotation rate f2 > can be taken to be that in which the sphere is at rest 
only if the assumption of infinite mass is made. For otherwise, the sphere will 
rotate at variable rate in response to the instantaneous amount of angular 
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momentum in the fluid component of tlie coupled system, to conserve the 
combined angular momentum. 

Recall that the vorticity is given by the Laplacian of the stream function, 

u = Atp 

In order to study the problem in a more tractable framework, we restate 
the problem in terms of local vorticity or lattice spins, following Lim [3]. 
The relative zonal and meridional flow is given by the gradient of the stream 
function and will be denoted as Ur and Vr respectively. In addition the flow 
due to rotation also contributes to the kinetic energy, it is entirely zonal and 
will be denoted as Up. The kinetic energy in the inertial frame will be the 
key objective functional in our work, 

H'[q] = dx[{Ur + Upf + v^] (1) 

= If dx[ul + vl + lUrUp + uj] (2) 

= -o / dxq^)^- / dxu\. (3) 

1 1 1 2 

The constant term ^ ||'Up||2 will be discarded below - we thus work with the 
pseudo-kinetic energy. 

Note that Stokes theorem implies that circulation on the sphere with 
respect to the rotating frame is zero, 

a; = (4) 

52 



4 The Discrete Model 

In the above discussion the relative vorticity is a function on S"^ given by 
the Laplacian of the stream function. [TT] , [12] , [3] , [1] utilize a transforma- 
tion to restate the coupled fluid-sphere model in terms of a relative vorticity 
distribution on the sphere which is then approximated over a discrete lat- 
tice. Proceeding from this transformation we consider a model defined by 
a uniform distribution of nodes {xj} on the sphere; edges spanning a pair 
of sites Xi.Xj are assigned interaction energy J(ij)- The kinetic energy of 
the fluid component is then given by pairwise coupling of local vorticity and 
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coupling of local vorticity to the planetary spin. Without loss of generality - 
for a single step renormalization based on block spins later relates this Ising 
model to a more physical model with real valued spins - we consider only 
local vorticity with spin of +1 or —1. As these spins are intended to model 
fluid flow, Stokes' theorem requires for 2A'" vortices of equal magnitude on 
the sphere that N vortices be of +1 spin and N vortices be of —1 spin. Under 
this restriction we consider the state space with the objective of finding the 
statistically preferred states. 

Formally, the state space consists of functions of the form 

uj{x) = ^Hi{x)si 

i 

where 

An 

Hi = - x). 

The vorticity Sj is assigned values ±1 which models the direction of the spin 
at a node. Thus we have the following simple form for the kinetic energy of 
the system, as a discrete approximation to the pseudo-kinetic energy, 

Hn = - ^ 'J{i,j)SiSj - ^ FiSi (5) 

where Y^' represents a summation with diagonal terms excluded, 

J(,,) = I f dxH,{x)G{H,){x) (6) 

Fi ^ [ dx cos e{x)G{Hi){x) (7) 

and G{-) is the inverse of the Laplace-Beltrami operator on the sphere. 

The usefulness of this form is clear; it closely resembles the lattice models 
of solid state physics, and becomes accessible to the rich theory of methods 
of that field. In the intensive or continuum limit we see that the coordination 
number of each vortex diverges. The model is not long-ranged in the technical 
sense, however, as the sum of interaction energies is finite. This differs from 
the Ising model in that every site is connected to every other site in the graph, 
the sum being over all lattice site pairs {x, x'). Note that as f2 ^ the spin 
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lattice model becomes the pseudo-kinetic energy for the non-rotating regime. 
It is easily seen that, 

•^(m) = (8) 
= -—cose{x,). (9) 

5 Bragg- Williams Approximation 

The central feature of the Bragg method is the approximation of internal 
energy of a state by its long-range order We will have to alter the previous 
Bragg method, as edges in the model have different energies as a function of 
their length, and furthermore, the rotation of the planet adds a kinetic energy 
term analogous to an inhomogeneous magnetic field. The method presented 
in this paper estimates internal energy from local order over domains on the 
sphere. The implicit assumption in this method, as in the original Bragg 
method is that the distribution of types of edges is homogeneous between 
domains. Specifically this is done by defining a partition of the sphere, into 
domains or blocks labeled {^} which contains in principle many original 
lattice sites. For each domain ^ we define notation N^{N^) = the number of 
sites in cell ^ which are positive(negative). Note = + A*"^. We define 
for every partition element ^ the local order parameter as: 

a, = 2^ - 1 (10) 

The method consists of approximation of important quantities by using - 
instead of the original discrete spin values - a local probability of spin value 
based on the local order parameter defined by the above block averaging 
procedure. Then the probability of any spin s in domain ^ being positive is 

P+ = Prob{s = +l\s e^} = ^-^, (11) 

and the probability of the spin being negative is 

p- = Prob{s = -l\se^} = ^^^^. (12) 

The energy is the coupling of spin domains, and entropy is calculated via the 
Shannon information entropy of spin variables {si}. 
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5.1 Statement of Equations 

The notion developed above of local order leads to a simple method of quan- 
tifying interaction types. Specifically, pairwise interaction is dependent on 
several parameters; the above interpretation lead to derivations of equations 
relevant to the renormalized spin-lattice model. Pairwise interaction occurs 

between all sites in three types {ss') = (++), ( ) and (H — ). Depending 

on the respective domains of the interacting sites, we use probabilities of spin 
distribution - coarse-grained variables - to inform probabilities of spin interac- 
tion. Below the probability P^r = Prob{si = +1, Sj = — G C,,Xj G , 
of an edge spanning domains ^ and ^' being of type (H — ), is calculated. 
Clearly, 

P++ + P^V + P^-f = 1 (13) 
which is identical to the relation 

n7 + - = 1 - 2^1" (14) 
The LHS of (fT^ arises below in (!20l) . and is the sole contribution of pair- 
wise order to the free energy. Thus we are content to calculate the order 
probability P^r only, as follows (■C' 7^ 0- 

Proh.., {+-} = = ^ ^ + — ^ ^ 



2 

Analogously we have, 



1 -|- cr^ 1 — 0"^/ ^ 1 + cr^/ 1 — (T^ 1 — cr^a^' 



Prob,, {+-} = ^ 

_ 1 + org 1 _ _ 1 - a| 



2 2 2 

for A^^ sufficiently large. 

Finally, we state the relation 

1 - 2P+r = a^a^' 

We have made this estimation based on the assumption that the number of 
spins in any domain ^ is sufficiently large, in fact in the thermodynamic limit 
of lattice models it is standard to take the number of spins to infinity. This 
must be done, however, on the finite surface of the sphere. 
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5.2 Estimation of Important Quantities 

In this section we estimate quantities in terms of the Bragg course-grained 
order for the system. However, first we will note that in the process of 
coarse-graining the spin values, enstrophy is now only bounded above, (t| < 
1. Although this diverges from other approaches where the enstrophy is 
constrained microcanonically, from our knowledge of the coupled barotropic 
fluid-sphere physical system - otherwise known as the Principle of Selective 
Decay for 2d flows - any bound on the relative enstrophy will be sufficient to 
make the resulting statistical mechanics model a well-defined one. A paper 
by the first author complementing this one uses a simple mean field 
method while enforcing enstrophy constraints in an averaged manner. 



Internal Energy The discrete approximation to the pseudo-energy (jS]) can 
be rewritten in the form, 

where edges are summed as members of domain pairs. The factor of | arises 
from double counting of domains. 

As consequences of the above assumptions - the distributions of positive 
and negative sites is homogeneous in a partition element - we calculate the 
Bragg internal energy in terms of the vector {cr^}. 



k,l \{i,j)\(xi,Xj)££,jX(,k / „ k \i\xie(,k 



where ()^ denotes Bragg averaging which will be made clear below. In view 
of the assumption of uniformity or homegeneity, we define the mean energy 
of an edge connecting two domains (or an edge from a domain to itself) in 
terms of area average 

Ks.i., = (ln(l-x-x')|(x,x') G^xO- (16) 

Further, define 

= - I ^cosQ{x) (17) 
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to be the area-weighted, mean couphng of a site with the external field. 
Returning to the energy functional, we focus on the individual terms, 

( Yl Jii,jy^i-^n (18) 

= ( J(i,j)H^~ s,Sj) - Y J{i,j)H^ + SiSj)) (19) 

\ii,j)\ixi,Xj)€ix^' i^,j)\(x^,XJ)e^x(' I ^ 

= (Pi' + p,v - { E %^)) (20) 

= a^a^,V^V^>K^^' (21) 

In ( |20|) we have the course-graining approximation over the edges in ^ x 
The expected interaction of a partition element with itself differs by a factor 
of I to account for double counting, 

J2 I (22) 

= \i Y J{i,j)5{'^- SiSj)- Y (1 + SiSj) y23) 

\(ij)|(xi,Xj)e5x5 {i,j)\(xi,Xj)i^£,x(, I ^ 

= + p,f - nf) { E (24) 

\{i,j)\{xi,Xj)(^(,x£, I 

= l^MK,, (25) 
Finally, the external interaction becomes 

FiSi \ = {si\xi E {Fi\xi G = j dx ^cos9(x) = cr^V^L^ 
Thus we have the simple expression for the Bragg internal energy 

u = -\Y. - ^ E ^^y^M^- (26) 
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Entropy It should be noted that the nature of our model is constrained to a 
constant sized sphere, and therefore no extensive quantities exist, only inten- 
sive ones. Since the coarse-grained or renormalized spins {cx^} are equivalent 
to the probability distributions in ffTTl) and f|T2l) . we therefore calculate -as 
is standard in equilibrium statistical mechanics- the Shannon entropy of the 
renormalized or coarse-grained state in terms of the local order parameters 

i 

This entropy is weighted by area of the domains. 



1 + (Tf . , 1 + (Tf. 1 

^ In ^' + - 



a.. , 1 



^5. 



(27) 



The Free Energy The Helmholtz free energy is then, 

^^J = U-TS 

C\ / J SI SI S^SJ / J SI Si 



1 + (Tf . , 1 + 0"f , 1 — CTt, , 1 — CTf . 

^' In ^ + ^ In ^ 



2 2 2 

constrained to the set of erg which are solutions to 



(28) 



The goal is to find critical points of the free energy with respect to {0"^^}^. 
Enforcing (128|) . we get a Lagrange multiplier problem. We must find the 
critical points, given by the simultaneous solution of the m equations. 



along with the constraint equation 



2 

Equivalently we write, = 1, .., m 



CTg. = tanh 



which is a m— dimensional fixed point equation for {ct^^}Y!=i- We will show 
that a well-defined continuum limit for this equation exists in the form a 
fixed point equation in Hilbert space L2(S'^). 
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6 Polar State Criteria 



In this section we implement the Bragg method on a two domain partition 
of the sphere. We find the order of the system varies continuously with the 
temperature and planetary spin. The system exhibits a continuous phase 
transition of second order at a negative critical temperature in the non- 
rotating case. The large spin f2 > thermodynamic regimes exhibit a con- 
tinuous phase transition at a positive critical temperature between a weakly 
counter-rotating ordered state and a mixed state. A contionuous transition 
at (3 = near-mixed counter-rotating states at high positive temperatures 
and near-mixed pro-rotating states at negative temperatures far from zero. A 
final continuous phase transition in the negative regime occurs when negative 
temperatures approaches zero to an ordered strongly pro-rotating state. 

The sphere is partitioned into the northern and southern hemispheres, 
labeled 1 and 2 respectively. In this case we have only parameters ai and a2 
representing local order parameters of the northern and southern hemispheres 
respectively. Then a2 = —cri from (jlj). This leads a simple approximation of 
the free energy 



:i + a,)lni±^ + (l-.,)ln^^' 

(29)" 



and fixed point equations for the extremal values 

ai = tanh [p{QLi + A) + (3aiVi{K^^i - K^^2)] 

-ai = tanh [-piQLi - A) - paiVi{Ki^i - ^1,2)] , 
which reduces easily to the fixed point equation in one variable 

ai = tanh [(3nLi + (3aiVi{Ki^i - fsTi.s)] . (30) 

Here Li = —1 and Vi = 2tt. A simple monte carlo integrator was used to 
estimate K = -V^ {K^^x - Ki^2) ~ 120. 



Non-rotating The non-rotating case is given in equation fl30|l by setting 
f2 = 0. For positive (3 the RHS of fl30|) is decreasing while the LHS increases, 
thus there exists a unique solution cxi = 0, which corresponds to the mixed 
state. 
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In the negative temperature domain /5 < we must maximize the free 
energy, according to an easy extension of Planck's theorem to negative tem- 
peratures. For this case it is clear that a solution is found at cxi = 0, cor- 
responding to a mixed state. Whether other fixed points exist depends on 
the slope of the RHS of fl30|) at ai = 0. Graphically it is clear this situation 
arises when the slope of the RHS has a slope of one, thus there is a critical 
quantity given by P'^Vi^Ku — K12) = 1. Thus for 

> /3 > = 

the point di = is the only stationary point (see figure ([1])). 




-1 -0.5 



0.5 1 



Figure 1: Left: The LHS and RHS for several values of (3 of equation (!30|) . 
Right: Corresponding plots of the free energy. For all plots = 0. 

For /5 < < there exist non-zero maximizers of the free energy cr~ < 0, 
(7~^ > 0; the tanh function is odd, thus cr~ = — cr~*'(figure ([T])). The free energy 
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is symmetric here so we need only consider fixed points cr"*" and cr = 0. 
Consider the Taylor series of the free energy around zero 

< $((t) - $(0) (31) 

= {K + 2Mb)(t^ + —TkB(T^ + 0{a^) (32) 

3 

We therefore have a polarized ordered state for /3 < /3° < 0, which corre- 
sponds to solid-body rotating flow. A standard symmetry-breaking phase 
transition is thus predicted by the fixed point equation in this Bragg mean 
field model, in agreement with both the simple mean field theory [19] and 
Lim's spherical model [12], [5], [I5], [1], [1]. 
We have established the result. 

Proposition 1: In the non-rotating case of the Bragg model, there is a 
negative temperature phase transition to a solid-body-rotating ordered state 
for j3 < (3^ = —2ti/K < 0. For all other values of the temperature - positive 
and negative - the most probable state is a mixed vorticity state. 



Rotating In the rotating regime the free energy can be approximated by 
a series as 

■qi^ ^ (a) = -AnTkB In 2 + 2^](T + {K + 2-KTkB) + ^^^a^ + O [a^] (33) 

3 

For > 0, and independent of whether the temperature is positive or 
negative, the RHS of the fixed point equation (|30ll is zero at 

-2tiVL , , 

-a, = < 0, (34) 

since K > 0. Moreover, it is easy to see that the zero ai satisfies 

- 1 < ai < (35) 

if and only if planetary spin Q - taken to be non-negative by convention - 
satisfies the condition 

0<Q<n,, = — (36) 
Ztc 

where Qcc > is independent of temperature - that is, planetary spins are 
not too large. This condition turns out to be significant for the determination 
of the fixed points of fl30l) in the rotating case which separates naturally into 
positive and negative temperature subcases considered next. 
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6.0.1 Positive temperatures 



In the case of positive temperatures, by virtue of the property that the RHS 
of fl30|) is decreasing in cji and the fact that it is bounded between ±1, there 
is always a fixed point o"i of fl30l) satisfying 

-1 < CTi < 0. 

By virtue of the shape of the graph of RHS of (!30|) . this fixed point also 
satisfies 

ai < ai < 0. (37) 
Furthermore, as /5 — > oo, cxi \ max ((Ti, —1). Since the zero ai satisfies 

- 1/2 < ai < (38) 
if the planetary spin Q is smaller than the threshold value 

n, = ^>0, (39) 

we deduce from equation (137|) that the fixed point also satisfies 

- 1/2 < ai < (40) 



for planetary spins smaller than Q^- From the fixed point equation fl30|) we 
find a phase transition where a < —1/2 at positive critical temperature 

4ttQ — K 
T+ (fi) kB = \ 
^ ' 27rln3 

Thus, because Qc is independent of the inverse temperature (3 > 0, we have 
established the following result where we have taken the cutoff |cri| = 1/2 to 
separate mixed states from ordered states, 



Proposition 2: (a) The most probable state in the Bragg model is 
the mixed vorticity state for all positive temperatures when planetary spins 
are smaller than - there are therefore no phase transitions in positive 
temperatures in this case. (This is shown in figure ^ where the fixed points 
are plotted vs T for very small Q < Q^-) (b) On the other hand for planetary 
spins that are not small, that is 

^>^c = ^, (41) 
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the fixed point G (—1,0) has continuously increasing long range order 

as temperature decreases (see figure\^ below) and helowT^iVt) > 0, the sta- 
tistical equilibria a\ < —1/2 is an organized counter-rotating physical flow. 

This is in agreement with the resuh obtained in the simple mean field 
approach in the case of positive temperatures [19]. Properties of the first 
threshold value of planetary spin are clearly shown in figure ([2]) and figure 3 
where fixed points and free energy respectively are plotted vs T for planetary 
spin Q = Qc- 
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Figure 2: A graph of a vs Tks when Q = fl^- 
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Figure 3: LEFT: A graph of the fixed point when 13 > 0,Q = Qc- RIGHT: 
The free energy 

6.0.2 Transition between positive and negative temperatures with 
large |r| 

For all values of planetary spins fl, the most probable state changes smoothly 
through mixed states between high positive temperatures T ^ 1 and large- 
absolute- valued negative temperatures T<^— 1. AtT^l, the preferred 
mixed state or fixed point has a slight negative angular momentum - it is 
slightly anti-rotating. For T ^ —1, it is the reverse - the fixed point is a 
mixed state with a slight positive angular momentum or pro-rotation bias. 
These facts are clearly shown in figures (jl]) and ([5]) for several distinct values 
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6.0.3 Negative temperatures 



We divide the analysis for negative temperatures into two natural categories - 
where the fixed point equation has the possibility of multiple fixed points and 
where it always have exactly one fixed point. In both categories, as negative 
T increases (or decreases in absolute value), and for all values of planetary 
spin, the slightly pro-rotating mixed state described in the last subsection 
gives way continuously to a strongly pro- rotating state at < - this 

is chosen to be the value of f3 for which the fixed point ai{(3~) = 0.5. This is 
clearly depicted in figures 2, 4 and 5 for several values of the planetary spin. 
Solving (1501) for Tks gives 

27rn + (tK , , 

T^^B = -TTCT- (42) 

Thus, 

However, figures 2, (j4]) and ([5]) also indicate that there is another tem- 
perature threshold when the spin is not too large, namely 

f3^ < (44) 

where < is the value for which multiple fixed points appear when {3 < 
and when VL < VLcc - both defined below. We can find an estimate for (5^. 
Critical value leading to the apearance of two more counter-rotating solutions 
of (!30l) is given by 

^^dTks^ K(l-a^)ln(i±g)+47r^] + a2ir 
da .(l-a^)(ln(l^))^ 

from which it follows 

An approximation of the numerator on the RHS of Il5]by a third order poly- 
nomial in a gives 



K 



(47) 



We can now conclude that 
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Lemma 3: The critical temperaturcT^ is increasing in Q 

Proof: Follows imeadiately from fHBl) and fB7|) 

One fixed point In the case of negative temperatures, the physical bound 
|cri| < 1, in the fixed point equation implies that there is another threshold 
on the planetary rotational rate. When 

n>n,, = ^^^^""^"^ = 2n,, (48) 

there are only pro-rotating solutions ai > to the fixed point equation 
because any negative fixed point af < - if they exist - must satisfy 

O-f < (Ti < 

where the zero ai must in turn satisfy 

(Ti < -1 

in view of fl35|) and fl36l) . which leads to the contradiction erf < —1. This is 
depicted in figure (jl]) for > Qcc- 
Thus, we have the result. 

Proposition 4: For large spins Q > Qcc, there is a single negative- 
temperature transition at [3~ {VL) < between the mixed and the strongly pro- 
rotating state for \T\ <^ 1. 

Multiple fixed points On the other hand, when Vl < Qcc, it is possible 
for the argument of tanh to be zero at values of ai within the physical range 
— 1 < 0"! < 0, making it in turn possible to have counter- rotating or mixed- 
state fixed points that satisfy 

-1 < CTi < CTi < 

if in addition, 

P<Pc< 0. 
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Figure 4: A graph of cr vs Tks when f2 > Vt^c 

In the case VL < VLcc-, and (3 < f3^, there are generically three fixed points 
due to the structure of tanh, one pro-rotating and in general two others - one 
of which is strongly counter-rotating and one mixed with a small counter- 
rotation - which merge into a single degenerate counter-rotating / mixed 
solution when /3 = This last threshold < is a consequence of the 
fact that —(3 determines the slope of tanh near and a large enough slope 
is required in order for the graph of tanh to intersect the line of the identity 
function. 

There is a transition at < between the mixed and the strongly 

pro-rotating state and another negative temperature transition < when 
Q < Qcc where aditional critical points arise in (l33ll . The Bragg free energy 
has a simple form in this case which we can exploit to easily see that the pro- 
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Figure 5: A graph of cr vs Tks when > ^ > ^■ 

rotating solution has greater free energy than the counter rotating solution 
in negative temperatures close to zero, as depicted in figure 6 for /3 > 

In the two domain case, > adds a linear term to the free energy 
functional, 

^T,n (a) = ^T,o (^t) + 2na (49) 

Recall from the nonrotating case that for > T > T° the free energy has 
two maximizers we will denote here as a~ and a^; the free energy in this 
case is even so = —a^. Therefore, from Lemma 3 and equation fH9l) . 

^T,n (^^) < ^T,n (fX") < ^T,f7 (^+) < (50) 

It is also of interest to note that < < < (T+ < o"+. In physical terms 
this implies that the ordered state > (in which positive relative vorticity 
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Figure 6: LEFT:Plots of the fixed points RIGHT: The free energy. Both 
plots are for the same value of /? < < 



dominates the northern hemisphere), is 'more ordered' than the symmetric 
solutions at = which is in turn more ordered than the counter-rotating 
fixed point erf < 0. We have demonstrated 



Proposition 5: For spins not too large, that is if Q < l^co then for 
hottest inverse temperature such that P < < 0, there are exactly three 
fixed points in the Bragg model but the most probable state is pro-rotating 
such that ai y 1. 
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6.1 Summary of main results 

In summary, the simple two domains case of the fixed point equation predicts 
that in the rotating problem, there are two critical values of the planetary 
spin, Qc in fl39l) and flee = '2flc- Below Qc, there is no transition at positive 
temperatures and the most probable barotropic fiow state is mainly a mixed 
state fHU]) with some small amount of counter-rotation which vanishes as 
temperature increases. For planetary spins above Qc, there is a continuous 
positive temperature transition from mixed states at T > T+(i7) to strongly 
counter-rotating barotropic states at lower positive temperatures. 

For all values of spin, there is a smooth transition through mixed states 
at /? = between slightly anti-rotating mixed states for T 3> 1 and slightly 
pro- rotating mixed states for negative T with \T\ ^ 1. Again for all values of 
spin Q, there is a continuous transition at l3~{Q) < between mixed states 
for cooler negative /5 > P~{fl) and strongly pro-rotating states for hotter 
f3 < f3~[Q). For large spins Q > Qcc, this is the only possible transition at 
negative temperatures. 

For intermediate to small values of spin, Q < Qcc, there is a negative 
threshold (3^ at which the possibility for multiple (meta-stable) thermody- 
namic equilibria first arises. We argued that when there are multiple fixed 
points, the pro-rotating branch has the largest free energy and thus contin- 
ues to be the thermodynamically stable macrostate. We conclude that for 
all spins, there is a single negative- temperature- transition at /3~ {fl) < to 
the pro-rotating state for the hottest negative temperatures - those with the 
smallest absolute values corresponding to the largest kinetic energy. 

The complete analysis of the transitions in the non-rotating case and in 
the rotating case agrees well with the simple mean field theory and Lim's 
spherical model. 

7 The infinite dimensional nonextensive limit 

As noted above, the model does not have an extensive limit from which 
we can derive intensive quantities. We have therefore bypassed the standard 
thermodynamic construction and directly built an non-extensive model. Con- 
vergence theorems support our continuum model. In this frame work we have 
had no need to include stipulations on the types configuration of domains. 
In light of this, we will consider a sequence of domain systems which are 
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constructed from Voronoi cells of a uniform mesh as in Lim and Nebus[ 
This affords us an advantage that we may assume the coarse grained domains 
to be spatially symmetric. We henceforth refer to the sequence of domain 
systems and their associated ensembles as a Bragg process. The limiting en- 
semble of a Bragg process is the space ^2(5"^) of functions a : 5*^ h-* [—1, 1], 
which has an expression of the free energy analogous to the discrete case, 



+Tk 



B 



1 

dx 



dx / dya{x)a{y)K{x,y) — Vl I dxo{x)L{x 



1 + aix) , 1 + aix) 1 
— In — + — 



aix) 1 — aix) 
In 



from which we can easily recover the fixed point equation. The free energy 
must be extremized over the subspace of functions constrained by Stokes' the- 
orem S = |(T : S"^ 1-^ [— 1,1]|0 = '^a;cr(x)}. Another physical constraint 
namely bounded relative enstrophy - the square of the L2 norm of the coarse- 
grained relative vorticity field a{x) - is implicit in the condition that |cr| < 1 
in the fixed point equation which can be derived in the continuous case sim- 
ilarly to the above discrete case, 



a[x) 



tanh 



[3{nL{x) + X) + P / dya{y)K{x,y) 
52 



(51) 



It is clear from the definition of K f[T6|) that we have recovered the inverse 
Laplacian G and we can rewrite fl5T]) as 



AV' = tanh [13 {nL{x) + X + ^p)] 



(52) 



where is the course grained stream function. From this form, remembering 
that L{x) is proportional to the first spherical harmonic and (A^) = we see 
that A = 0. Indeed, for any ip we have A = implies (tanh [P {QL{x) + A + ip)]) 
and, 



d 

dX 



(tanh [j3 {nL{x) + X + ip)]) 



d 

dX 



tanh [p {nL{x) +X + ij)] 



= (/3cosh-2 [P {VlL{x) + a + ^)]> 

serves to show the average is monotonic in A. 

While the original discrete model was limited in its degrees of freedom, we 
now have a continuum of freedom for the spin since now the coarse-grained 
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state is given by a bounded function cr{x) in ^2(6'^) with zero circulation 
Jg2 dx(j{x) = 0. Further the Bragg estimation of interaction energy became 
very accurate in the continuum modeL The price paid, hke in all mean field 
models, is the artifact of the entropy estimation. However, the Shannon 
entropy of a binary variable is exactly what enables our derivation of the 
tanh fixed point equation above. The figure below shows the spin states of 
some thermodynamic regimes obtained numerically. 
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Figure 7: Some fixed point solutions 



8 Conclusion 

In conclusion it must be emphasized that the combined angular momentum 
of the fluid and solid sphere is conserved even when that of the fluid alone 
is not. Similarly, without loss of generality for the purposes here, we can 
assume that the complex torque mechanism responsible for coupling fluid 
to sphere is conservative, and thus, the combined energy of the fluid-sphere 
system is conserved. 

In taking the mass of the rotating sphere to be infinite here and in [12], [5], 
we have essentially modeled the statistical mechanics of the coupled fluid- 
sphere problem as one having separate reservoirs for angular momentum 
and energy. It was shown in [12] and [19] that the particular form of the 
kinetic energy functional for coupled barotropic flow on a sphere allows us to 
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conveniently combine these two reservoirs into one so there is only one inverse 
temperature (3. Introducing a distinct reservoir for angular momentum, that 
is, imposing a separate canonical constraint on the angular momentum of the 
fluid will not change the physics but will add to the clutter in the resulting 
model. 

We remark here that a series of experiments on decaying and forced- 
damped (or stationary) quasi-2D turbulence in square domains with no-slip 
and stressless boundaries - see and references therein - have highlighted 
the role of variable angular momentum in inverse cascades to self-organized 
structures. In particular, normal forces at the boundary of the square can 
in principle produce a torque without requiring shear or the presence of fric- 
tion at the boundary. This conservative mechanism for angular momentum 
transfer between fluid and boundary - despite the obvious mismatch in ge- 
ometry - provides the best experimental comparison that we know to the 
(unresolved) torque mechanism introduced and developed by the first author 
into the coupled fluid-sphere statistical mechanics models in a series of recent 
papers. 

For a more detailed view of the complex torque mechanism that couples 
the fluid to the rotating sphere, one will have to invoke the observed and 
theorized fact that a baroclinic component of the fluid forms Hadley cells 
driven by non-uniform solar insolation which in turn serves to transfer high 
momentum fluid in the equatorial region up to the middle atmosphere and 
mid-latitudes. The planetary boundary layer is a very complex region that is 
responsible for transferring the solid sphere's angular momentum to the fluid 
as it flows towards the equator from mid-latitudes and the fluid's angular 
momentum to the sphere between the polar region and mid-latitudes. Any 
attempt to model this torque mechanism explicitly will result in technical 
difficulties that can only be resolved by a huge computational effort. Fortu- 
nately, the very essence of statistical mechanics - namely coarse-graining - 
allows us to formulate a viable theory without explicitly resolving this com- 
plex torque mechanism. 

The results derived by using Bragg's method are in good agreement 
with those found through Monte - Carlo simulations of the spherical energy- 
enstrophy model by Ding and Lim [Ij, [4J and with a simple mean field theory 
[19j, in the non-rotating, rotating positive temperature and possibly rotating 
negative temperature problems. More details can be found in the recent book 



[TT] and on the website http:/ /www.rpi.edu/~limc[ Further work should be 
done, however, to refine the mult i- dimensional Bragg method and extend it 
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to more complex geophysical flows such as shallow water flows coupled to 
a rotating sphere. The inflnite-dimensional flxed point formulation in the 
previous section gave us an interesting nonlinear elliptic equation fl52|) that 
should be further analyzed. Complementing the work here, the flrst author 
recently presented exact solutions to the spherical model for the energy- 
relative enstrophy theory of the coupled barotropic fluid-sphere system [12], 
[5]. 
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